Evolution of collective N atom states in single photon superradiance 
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We present analytical solutions for the evolution of collective states of N atoms. On the one hand 
is a (timed) Dicke state prepared by absorption of a single photon and exhibiting superradiant decay. 
This is in strong contrast to evolution of a symmetric Dicke state which is trapped for large atomic 
clouds. We show that virtual processes yield a small effect on the evolution of the rapidly decaying 
timed Dicke state. However, they change the long time dynamics from exponential decay into a 
power-law behavior which can be observed experimentally. For trapped states virtual processes are 
much more important and provide new decay channels resulting in a slow decay of the otherwise 
trapped state. 



The long standing problem of collective spontaneous 
emission from N atoms 111, Jsjl 1S a subject of growing 
recent attention USIEIEHIaQl and debate [nj. Effects 
of virtual processes are of particular current interest @, 

EMM- 

Here we consider a system of two level (a and b) atoms, 
E a — Eb = fko. Initially there are no photons and one of 
the atoms is in the excited state a, but we don't know 
which one. That is the system is prepared in a collec- 
tive N— atom state. The question then is how such a 
collective state evolves with time. 

Atoms interact with common electromagnetic field and 
the interaction Hamiltonian is given by 



Hi nt = 



N 

k j=i 



+ a k e 
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where Oj is the lowering operator for atom j, Ok is the 
operator of photon with wave vector k, is the atom- 
photon coupling constant and Yj is the radius vector of 
the atom j. Evolution of the atomic system is described 
by the state vector 
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where |&i &2 -^iv > is a Fock state in which atom j 
is in the excited state a and all other atoms being in 
the ground state b. We disregard polarization effects, 
that is treat photons as scalar and assume that initial 
state evolves slowly compared to the time of photon flight 
through the atomic cloud (the opposite limit has been 
studied in Q). 

Decay of an initial state occurs via real and virtual 
processes in which a virtual photon is emitted and then 
reabsorbed. In particular, due to counter-rotating terms 
in Hamiltonian {T]) virtual processes couple the single- 
atom excited states with those in which two atoms are 
excited. If all virtual processes are taken into account 



then for a dense atomic cloud evolution of the system 
is described by an integral equation with an exponential 
kernel @, || 



d(3(t,r) . N f ,exp(ifc |r-r'|) 



Of 
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-my), (3) 



where V = 4ttR 3 /3 is the volume of the spherical atomic 
cloud, fc — uj/c and 7 is the single atom decay rate. We 
assume that atoms are uniformly distributed with density 
N/V in a sphere of radius R. 

If we ignore virtual contributions then Eq. ([3]) reduces 
to an equation with sinusoidal kernel 



d/3(t, r) = _ N r , sin(fc |r- 
dt 1 V J r fco|r-r 



r'|) 



my)- (4) 



Here we solve Eqs. ^ and (U]) analytically for two initial 
conditions, namely the |+ > "timed" Dicke state 



= e ik °- r , 



(5) 



which is prepared by absorption of a single photon with 
wave vector ko (fco = w/c) [3, 0|) and the symmetric 
Dicke state [l[ 



/3(0,r) = l. 



(6) 



For a large atomic sample R 3> A (A = 2irc/ui is the 
wavelength of the emitted photon) the |+ > state ([5]) is 
superradiant, while ([6]) is a trapped state undergoing very 
slow decay. As we show below, virtual processes yield a 
small (yet interesting) effect on evolution of the rapidly 
decaying |+ > state. Such states decay mainly via real 
Weisskopf-Wigner spontaneous emission processes. How- 
ever, virtual processes can substantially modify the dy- 
namics of trapped states and provide a main channel of 
decay. 

Figs. [T]|4] summarize our main findings. For a small 
atomic cloud R -C A symmetric state © exponentially 
decays according to Eq. ((4| with rate T = ^7 without 
coupling to other states. This result has been obtained 
by Dicke Our Figs. [1] and [2] show, however, that 
virtual processes excite other states with a few % prob- 
ability even in the small sample (Dicke) limit. For a 
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large cloud Eq. ((4]) predicts that symmetric state ((6|) is 
trapped, however, virtual processes lead to its slow decay 
as shown in Fig. [3j On the other hand, evolution of the 
rapidly decaying |+ > state is only slightly affected by 
virtual processes (see Figs. Q]and|4|). For a large sample 
such processes excite other states with less then about 
10% probability. Thus the timed Dicke state ([5]) is, to a 
good approximation, described by Eq. (f5|) which ignores 
virtual transitions. However the symmetric state J6]) is 
strongly effected by virtual processes as per Fig. [3] 
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FIG. 2: Probability that atoms are in the symmetric state 
P s {t) obtained using the exp (solid line) and sin (dash line) 
kernels for R = 0.01 A. Initially atoms are in the symmetric 
state © and T — Nj. Insert shows probability to find atoms 
in any other state but symmetric state calculated for R = 
0.01A (solid line) and R — 0.03A (dash-dot line) from equation 
with exp kernel. 



FIG. 1: Probability that atoms are excited P(t) for small and 
large (insert) atomic clouds calculated using the exp (solid 
line) and sin (dash line) kernels. Initially atoms are prepared 
in the |+ > state, F = N-y (for U«A) and F — 3N-y/2(k R) 2 
(for R ^> A). For small atomic sample the |+ > state and the 
symmetric state (|6]) are the same. 

Next we discuss the evolution of |+ > state in detail. 
For R^> X equation with sin kernel ((3]) gives 0, 0] 
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where 



r 



3Nj 



(7) 
(8) 
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Here we hnd that equation with exp kernel ([3]) yields 

(3(t,r) =e* k °- r lf(t 7 r)+ig(t,r)cos6}, (9) 
where 9 is the angle between ko and r, 



f(t,r) 



9(t,r) 
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FIG. 3: Probability that atoms are excited P{t) obtained 
using the exp (solid line) and sin (dash line) kernels. Initially 
atoms are in the symmetric state (|6]), R = 5 A and Y is given 
byEq. ©. 



and Jo (2) is the Bessel function. Answer (|9|) is universal 
in the sense that state evolution is determined by the di- 
mensionless time Tt and change of the sample size simply 
results in time rescaling. 
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Next we calculate the probability that atoms are ex- 
cited as a function of time 




P(t) 



1 

V 



dr\(3(t,r)f 



For the integral equation with sin kernel 



(17) 



(18) 



For P(t, r) given by Eq. ([9]) one can calculate the integral 
in Eq. (117[) numerically for any t, while for f > 1/T we 
find 



7V2 



0.21 



ISTrVrt 



Ft 



(19) 



FIG. 4: Probability that atoms are in the |+ > state P+(t) 
obtained in the large sample limit from equation with exp 
(solid line) and sin (dash line) kernels. Insert shows proba- 
bility that atoms are in any other state but |+ >. Initially 
atoms are in the |+ > state and F is given by Eq. ©. 



To obtain solution ^ we used the identities 
exp(ifc | r — r '|) 



fc |r - r'| 



^(2m + l)P ro (f-r')x 



m=0 



Insert of Fig. Q] shows P(t) obtained using the exp ker- 
nel (solid line) and Eq. (TS]) (dash line). At t < 1/T 
the function P eX p(t) decays as e _2rt , while for t > 1/T it 
becomes closer to its asymptotic expression (fT9l . During 
the major part of the decay curve P eX p(t) exhibits expo- 
nential behavior (|18p and, thus, virtual processes have 
essentially no effect. However, virtual processes modify 
Pexp(i) at large time yielding the power-law decay (fT^jl . 
Such an interesting, although small, effect can be ob- 
served experimentally. 

For solution © the probability that atoms are in the 
|+ > state is given by 



Jm 

(k Q r')h ( A\k r), r>r' 

J in 

(k r)h ( A } (k r'), r<r> ' 



(12) 
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[A~6Tt+{Tt) 2 ]J 1 {2V2Tt) 



exp(ik • r) = ^2 *™( 2n + l)jn(k r)Pn(ko ■ f), (13) which for t > 1/T yields 



4n 
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dfl r ,P m (f-r')P n (k -f') = 6 mn ——P n {k -f), (14) 

In + 1 



where f and fco are unit vectors in the directions of r 
and ko respectively, P n are the Legendre polynomials 
and jn{ z ), hn\z) are the spherical Bessel functions. In 
the large sample limit the ansatz © yields the following 
equations for the slowly varying functions / and g 



4n{Tt) 3 / 2 



cos 



(20) 



(21) 



df(t,r) 
dt 



dr'f(t,r')-^ I <Ir' 1/(1 . ). (15) 



^M = £ / dr >f { ty), 

dt R 1 J v ' ; ' 



(16) 



with the initial conditions f(0,r) — 1 and g(0,r) = 0. 
One can solve Eqs. (fT5|) and (jT5J) using the method of 
Laplace transform which yields the answer (|10p and (fTTj) . 



In Fig. 0]we plot P+ p (i) obtained from Eq. (J2SD (solid 
line) and compare it with those found from equation with 
sin kernel P^ n {t) = e~ 2rt (dash line). The two curves are 
very close to each other. This means that virtual pro- 
cesses practically do not change evolution of |+ > state 
if it is considered separately. Without virtual processes 
the |+ > state directly decays into the ground state by 
emitting a photon. Virtual processes yield an extra de- 
cay channel in which energy is partially transferred into 
other atomic states. However, as one can see from Fig. 
HI the net decay rate of the |+ > state into all chan- 
nels remains practically the same with or without virtual 
processes. Insert shows probability that atoms are in any 
other state but |+ >. This curve demonstrates that dur- 
ing the system evolution the other states are excited with 
probability less then about 10% and, therefore, the effect 
of virtual processes is quite small for fast decaying states. 
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Next we discuss evolution of the symmetric state ([6]) . 
For such initial condition Eq. (U) with sin kernel can be 
solved analytically for any size of the atomic sample and 
yields 



sample limit 



where 



and 



(3(t,r) 



F = 



for l j 

koRcos(koR) — sm(foR) 
foR. — sin(feo-R) cos(fc i?) 



r = 



3 7 iV 
2ifoR) 2 



1 - 



sm(2k a R) 



2k R 



(22) 



(23) 



(24) 



Eq. (|22j) shows that at the beginning the atomic system 
decays with the superradiant rate (|24[) but quickly ends 
up in a trapped state 



0(r) = 1 + 2F 



sm(kor) 
for 



(25) 



Function (|25p vanishes in the small sample limit foR <C 
1, however, for large sample /3(r) ss 1 and state ([6]) is 
completely trapped. Probability that atoms are excited 
is given by 



6 [k Rcos{k Q R) - sm(k R)} 2 [l 
foR — sin(/coi?) cos(foR) (foR) s 



-2rti 



Pit) = 1 - 

(26) 

For a large atomic cloud R 3> A the evolution Eq. ([3]) 
with initial condition © can be also solved analytically 
and the answer is expressed in terms of the Bessel func- 
tions. In Fig. [3] we plot probability that atoms are ex- 
cited P(t) obtained from equation with exp (solid line) 
and sin (dash line) kernels. Initially atoms are prepared 
in the state ©■ Size of the atomic sample is R = 5 A. 
Insert shows behavior of P(t) for exp kernel on a large 
time scale which exhibits interesting plateaus and oscilla- 
tions. For t less then a few 1/F two curves are identical. 
For such time the real processes dominate and the initial 
state evolves into the state (l2"5j) which is trapped if we 
omit virtual processes. Virtual processes, however, re- 
sult in state decay as shown by the solid curve. State © 
overlaps with many eigenstates of Eq. Q. Eigen- 

states which decay faster contribute to evolution at small 
time. As time increases P(t) decays more slowly. How- 
ever, eigenfunctions of Eq. ^ are not orthogonal and, 
in addition, have different collective Lamb shifts. This 
makes state evolution richer. 

In the small sample limit R A the initial states ([5]) 
and © are the same. Equation with sin kernel gives 
Dickc result 



(27) 



where T = Nj. For equation with exp kernel the state 
evolution can be obtained by noting that in the small 
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are eigenfunctions of Eq. ([3]) with eigenvalues Q 



A„ 



96Nj 



n 2 (2n+l) 2 foR 7r 4 (2n+ l) 4 
Using the identity 

(-1)" 



n = 0,1,2, 



(29) 



1 = -E 

rrrrp f J 



irx ^ (2n + l) 2 

n— v ' 



sin[(2n + l)a 



one can expand the initial condition /3(0, r) = 1 in terms 
of (3 n (0,r). As a result, time evolution of the symmetric 
state is given by 
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R_ 



n . , 8R ^ (-1)" . 
P(t,r ) = — > 7^ — -TToSm 
v ' 7T T ^ (2n + l) 2 

n— / 



and probability to find atoms excited is 



P(A - ££ exp [-2Re(A n )t] 
nj ~7r 4 ^ (2n+l) 4 ' 



-A„t 



(30) 



(31) 



Fig. Q] shows Pit) given by Eq. (13~1]) (solid line) and 
compares it with the answer obtained omitting virtual 
processes P(t) = exp(— 2Tt) (dash line). The two curves 
are close to each other, but Eq. (|3"T|) yields a few % of 
population trapped which slowly decays with time. 

In Fig. [2] we plot probability that atoms are in the 
symmetric sate ^ obtained for R — 0.01A from Eq. (|3"0|) 
(solid line) and compare it with P(t) = exp(— 2Tt) (dash 
line). The two curves are very close meaning that the 
net decay rate of the symmetric state into all channels is 
the same with or without virtual processes. Insert shows 
probability to find atoms in any other state but symmet- 
ric sate © for R = 0.01A (solid line) and R = 0.03A 
(dash-dot line) obtained from Eq. (l30|) . Dependence of 
the imaginary part of A„ (collective Lamb shift) on n is 
the reason for oscillations. Period of oscillations is pro- 
portional to foR. The other states are excited with a 
few % probability. Thus, in the small sample limit, vir- 
tual photons also yield a small (but interesting) effect on 
evolution of fast decaying states. 

In summary, we consider evolution of two collective 
states of N atoms, the |+ > state which decays fast and 
the symmetric state which is trapped for R 3> A. We 
obtain analytical formulas for the atomic state vector as 
a function of time. We show that virtual processes yield 
a small effect on evolution of the rapidly decaying states. 
However, they change the long time dynamics from expo- 
nential decay into power-law which can be observed ex- 
perimentally. For trapped states virtual processes qual- 
itatively modify state evolution. Namely, they provide 
new decay channels which ultimately result in a slow de- 
cay of the otherwise trapped state. 
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